Thermal shifts and intermittent linear response of aging systems 
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At time t after an initial quench, an aging system responds to a perturbation turned on at time 
f w < t in a way mainly depending on the number of intermittent energy fluctuations, so-called 
quakes, which fall within the observation interval (i w ,t] [Sibani et al. Phys. Rev. B, 74, 224407 and 
Eur. J. of Physics B, 58,483-491, 2007]. The temporal distribution of the quakes implies a functional 
dependence of the average response on the ratio t/t w . Further insight is obtained imposing small 
temperature steps, so-called T-shirts. The average response as a function of t/t w ,eff, where t Wie g is the 
effective age, is similar to the response of a system aged isothermally at the final temperature. Using 
an Ising model with plaquette interactions, the applicability of analytic formulae for the average 
isothermal magnetization is confirmed. The T-shifted aging behavior of the model is described using 
effective ages. Large positive shifts nearly reset the effective age. Negative T-shifts offer a more 
detailed probe of the dynamics . Assuming the marginal stability of the 'current' attractor against 
thermal noise fluctuations, the scaling form £ w ,eff = iw, and the dependence of the exponent x on 
the aging temperatures before and after the shift are theoretically available. The predicted form of 
x has no adjustable parameters. Both the algebraic scaling of the effective age and the form of the 
exponent agree with the data. The simulations thus confirm the crucial role of marginal stability in 
glassy relaxation. 

PACS numbers: 65.60.+a, 05.40.-a, 61.43.Fs 
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I. INTRODUCTION 

Rejuvenation, memory and other intriguing aspects of 
out-of-equilibrium thermal relaxation, a process widely 
known as aging, continue to attract experimental and 
theoretical interest 0, 0, i, S H, S, 0, S S M HO, M 

m Hi HE IE HI- 

One starting point for theoretical analyses is the real 
space morphology of spin-glasses with short-range inter- 
actions, where thermalized domains form, with a linear 
size growing in time as R(t) oc t A @, 0, @, EH • The do- 
main size corresponds to the thermal correlation length 
and can be extracted from numerical simulation data. 
Furthermore, the temperature dependence of the expo- 
nent A can be rationalized by assuming that the energy 
barrier controlling domain growth increases logarithmi- 
cally in time [i, [la ] . The same logarithmic time depen- 
dence characterizes hierarchical models, which account 
for e.g. memory effects using configuration space, or 
'landscape' properties as starting point 0, [H, |2(|. A 
third theoretical approach focuses on the relationship 
between linear response and autocorrelation functions. 
Numerical [U H2] and experimental [H, H3, Hj| results 
show that the Fluctuation Dissipation Theorem (FDT) 
is applicable for observation times i b s < t WJ where t w 
is the time at which the perturbation is applied. While 
the FDT is broken for larger i b s , the proportionality be- 
tween conjugate response and correlation functions may 
be restored asymptotically for large t w and i b s , whence 
a general description flows in terms of effective tem- 
peratures [H, Hy, H3, HI]. Last but not least, fluctu- 
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ation spectra [22|, |29|, l30|, |3lL |32|, |33j indicate that ther- 
mally activated aging processes, linear response functions 
included, are controlled by intermittent transitions be- 
tween meta-stable attractors which irreversibly release 
the excess energy trapped in the initial configuration. 
These events, dubbed quakes, have statistical properties 
suggesting that they might be tri gger ed by thermal fluc- 
tuations of record magnitude |34l |35| . 

Based on the idea that metastable attractors of glassy 
systems are marginally stable, and that the transitions 
between them are hence induced by record sized ther- 
mal fluctuations, a 'record dynamics' scenario (3o| pro- 
vides approximate analytical formulae for several physi- 
cal quantities. The derivations ignore quasi-equilibrium 
properties, but are sufficiently general to support the 
universal character of off-equilibrium aging phenomena 
H HI H [H, EH EH. To buttress a comprehen- 
sive description of non-equilibrium aging based on record 
dynamics, this paper discusses numerical simulations of 
an Ising model with plaquette interactions, which has a 
glassy dynamical regime at low temperatures. The work 
extends numerical investigations (3ll . [32j] of the same 
model, and experimental and numerical studies of linear 
response in spin-glasses [22, 30]. 

Aging under isothermal conditions is considered first, 
and analytic formulae for the average magnetization [22l . 
l30l ] are verified. Secondly, the effect of 'T-steps' or T- 
shifts', small temperature steps of either sign applied con- 
currently with the magnetic field is analyzed in terms of 
effective ages. Replacing the age t w as a scaling param- 
eter, the effective age makes the response look similar to 
the isothermal response at the final temperature. Effec- 
tive ages give an excellent parameterization of the av- 
erage response for negative shifts, and a reasonable one 
for positive shifts. Fluctuation spectra cannot be fully 
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described by an effective age. 

For negative T-shifts, a n alg ebraic relation between ef- 
fective and true age exists [17| , and the analytical depen- 
dence of the corresponding exponent on the temperatures 
before and after the step has no adjustable parameters. 
Both the algebraic scaling and the for of the exponent are 
confirmed by the simulations. The final discussion em- 
phasizes the ramified connections of the present approach 
to other descriptions of aging dynamics. 



II. MODEL AND METHOD 

In spite of its simplicity and, in particular, in spite 
of its lack of quenched disorder and the triviality of 
its ground states, the p-spin model has full fledged ag- 
ing behavior (3li l42l . |43| . In the model, N Ising spins, 
Ui = ±1 are placed on a cubic lattice with periodic 
boundary conditions. They interact through the plaque- 
tte Hamiltonian 



H 



<Ti<Tj(Jk(Ji + Hi](t w - i)2jeri, (1) 
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where first sum runs over all the elementary plaquettes 
of the lattice, including for each the product of the four 
spins located at its corners. The second sum describes 
the coupling of the total magnetization a i to an exter- 
nal magnetic field. As expressed by the Heaviside step 
function r/(t w — t), the field jumps at t = t w from zero to 
H. 

The simulations utilize the Waiting Time Algorithm 
(WTM) a rejectionless algorithm endowed with an 
'intrinsic' time unit approximately corresponding to one 
Monte Carlo sweep. By choosing a high energy random 
configuration as initial state for low temperature isother- 
mal runs, an instantaneous thermal quench is effectively 
achieved. All simulations temperatures are within the 
model aging regime, 0.5 < T < 2.5. To improve the 
statistics, PDF data are collected over 2000 independent 
runs for each set of physical parameters. Other statistical 
data are collected over 1000 independent runs. 

In the following, the symbol t stands for the time 
elapsed from the initial quench (and from the beginning 
of the simulations). The field is switched on at time t w , 
and i bs = t — t w is the 'observation' time, during which 
data are collected. The external field jumps from zero to 
H = 0.3 at t = i w . The thermal energy is denoted by E, 
the magnetization by M and the fluctuations by SE and 
SM respectively. The latter are calculated as finite time 
differences over time intervals of length St <C £ bs- The 
average magnetization per spin is denoted by /kzfcm- 



III. RESULTS 

That intermittent physical changes occurring in an ag- 
ing process should be statistically subordinated to the 



quakes [22], |3CJ, [4CJ, |45j is a crucial element of a record 
dynamics description. For completeness, we therefore 
include recent supporting evidence [13] related to the 
present model. 

The circles in Fig. ([T]) (see Ref. [Hj]) show, on a loga- 
rithmic vertical scale, PDFs of the magnetic fluctuations 
SM in the interval [t w ,t w + £ bs]- The data in the left 
panel are spontaneous magnetic fluctuations (no field), 
while those in the right panel are effected by an external 
magnetic field switched on at t = t w = 1000. In both 
panels, intermittent wings extend a zero-centered cen- 
tral Gaussian spectrum. The magnetic field enhances the 
positive wing, reduces the negative one, and leaves the 
Gaussian spectrum unchanged. In the right panel, a con- 
ditional PDF (stars) is also shown. The PDF is obtained 
by removing from the statistics all magnetic fluctuations 
occurring within the same or the next St as energy fluc- 
tuations of magnitude SE < — 5. As the threshold cho- 
sen delimits the intermittent behavior of the heat flow 
PDF [32j , the filtering removes the magnetic fluctuations 
concurrent with the quakes. Since the resulting condi- 
tional PDF is nearly Gaussian, the intermittent magnetic 
fluctuations carrying the average magnetic response are 
all closely associated to the quakes. Furthermore, the en- 
ergy, and, in particular, the quakes, are hardly affected 
by the magnetic field. Therefore, the quakes dissipate the 
excess ener gy e ntrapped in the initial (zero field) config- 
uration [3l], 32] , and the field probes but does not modify 
the spontaneous intermittent dynamics. 

As argued in ref. [361 ] . the number of quakes in the 
interval (t, t') is theoretically described by a Poisson dis- 
tribution with average 



ni(t',t) = a{N) \n(t/t'). 



(2) 



In the theory, a(N) depends linearly on the system size 
N and not at all on temperature. The present model has 
nearly the same behavior [3^] • 



A. Isothermal aging 

For physical processes subordinated to statistically in- 
dependent quakes, the number, n, of quakes falling in a 
given observation interval assumes the role of the time 
variable in a homogeneous Markov process. The corre- 
sponding moments represent physical observables, and 
admit eigenvalue expansions. By averaging such expan- 
sions over the Poisson distribution of n, analytic formulas 
are obtained for the moments [12, H(J 5H which depend 
on the ratio t/t w between the end-points, t w and t, of 
the observation interval. These expressions can usually 
be truncated after a few leading terms. E.g. the average 
linear response of the p-spin model is well approximated 
by 



Mzfcm — b + a r 



ln(f ) + Mr ) A " 



(3) 
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FIG. 1: (color on line) (a): The PDF of the spontaneous magnetic fluctuations has a Gaussian central part and two symmetric 
intermittent wings. Data are sampled in the interval t w ,t w + t b s . (b): The same quantity (outer graph) when a field H = 0.3 
is turned on at i w = 1000. The left intermittent wing is clearly reduced, and the right intermittent wing is amplified. The 
inner, almost Gaussian shaped, PDF is the conditional PDF obtained by excluding the magnetic fluctuations which happen 
in unison with the quakes. The lines are fits to a zero centered Gaussian and to two exponentials, which are obtained using 
different subsets of the data. 



where bo, a m , b m and A m are parameters. The exponent 
X m is negative, whence, asymptotically, the behavior be- 
comes logarithmic. Figure ^ shows the average linear 
response versus t/t w for the temperatures indicated in 
the four plots. Each plot displays data sets (dots) for 
t w = 100,200,400,600,800, 1000 and 2000. The full line 
is a fit to Eq. © . In order to include the stationary con- 
tribution to the magnetization, which is ignored in the 
theory, the constant bo in Eq.© is modified by an ad- 
ditive £ w dependent term. The correction increases with 
t w , but stays within a few percent of the total. The black 
line is obtained by fitting to Eq. ([3]) . For the smallest val- 
ues of the abscissa, the quality of the data collapse, and 
consequently, of the fit, is higher the higher the tempera- 
ture. For larger values of the abscissa, all fits are equally 
satisfactory. In summary, the isothermal aging behavior 
is rather well described by Eq. (J3j). 



B. The effect of T-shifts 

The response function of spin-glasses [j| [46[ is sys- 
tematically affected by a small temperature change, a so 
called T-shift, applied during the aging process, usually 
together with a field switch at age t w . E.g. the logarith- 
mic derivative S(t w ,t b s ) of the average magnetization 
as a function of i bs peaks at i bs = t w under isother- 
mal conditions, and the peak moves to a lower, respec- 
tively higher, value when positive, respectively negative, 
temperature shifts are applied. As the same effect can 
be achieved isothermally at the final temperature using 
a shorter, respectively longer t w , an effective age t w>e s 
expectedly provides a natural and general parameteri- 



zation. A study of the intermittent heat transfer of the 
E-A spin glass model [T3| under a T-shift, and the present 
study, confirm and to some extent qualify this conclusion. 

Below, we consider how both negative and positive T- 
shifts affect the average response function of the plaquet- 
te model. We furthermore discuss how negative shifts 
affect the spectrum of intermittent magnetization fluctu- 
ations. Figure ([3]) shows the average magnetic response, 
vertically adjusted by a small t w dependent term as pre- 
viously discussed, for a negative T-shift imposed con- 
currently with the field switch at t w . Each panel of 
Fig. © contains three curves. In the lowest curve the 
T-shifted magnetization is plotted versus the scaling vari- 
able t/t w — 1. The two, nearly overlapping, and higher 
lying, curves are i) the T-shifted data again, now plotted 
versus t/t WiC g — l, with £ w , c ff chosen to maximize the data 
collapse, ii) Isothermal magnetization at the final tem- 
perature, plotted versus t/t w — 1. Clearly, scaling by a 
suitable effective age accounts rather well for the average 
linear response following a negative T-shift. Figure Q 
is similar to Fig. ([3]), except that the applied T shifts 
are now positive. In each panel, the highest lying curve 
shows the T-shifted magnetization versus t/t w — 1. The 
two lower lying curves are i) the T shifted magnetiza- 
tion data, now plotted versus i/i Wje ff — 1, and, ii) the 
magnetization curve obtained for isothermal aging at the 
final temperature, plotted versus t/t w — 1. Note that the 
transient effects produced by a positive shift are not fully 
described by an effective age. Consequently, the 'best' 
value of £ w ,efl is obtained by a fit which only maximizes 
data overlap for large t/i w values. 

Figure ([3]) shows that positive T shifts lead, as ex- 
pected, to effective ages smaller than the actual ages. 
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FIG. 2: (Color on line) For each temperature indicated, the four panels depict scaling plots of the average magnetization per 
spin, vertically adjusted by the constant &(i w ). The values of t range from t = i w to 400i w , with i w = 50,100,200,500 and 
1000, (green, blue, magenta, black, and red circles respectively). The black line is a fit to Eq.©- 



The best achievable data collapse is not as good as for 
negative shifts, especially for small values of t/t w . Since 
a positive T-shift destroys, totally or in part, the current 
configuration, it apparently induces additional dynami- 
cal effects not well represented by an effective age. Qual- 
itatively, the numerical results for positive shifts concur 
with a crucial hypothesis of record dynamics, namely 
that the attractors successively visited are all marginally 
stable, i.e. they are destabilized by a record sized thermal 
fluctuation. Increasing the temperature leads to larger 
thermal fluctuations, and hence to the destruction of the 
configuration reached at t w . Negative shifts have more 
subtle consequences, which yield to an analytical descrip- 
tion. We first consider how negative shifts affect the 
spectrum of the magnetization fluctuations. Secondly, 



the algebraic relationship between true and effective age 
is explained theoretically and checked against simulation 
data. 

A system undergoing a negative T-shift T = 1.9 — » 
T = 1.5 at t w = 1000 has, according to Fig. (O, second 
panel, an effective age i WjC ff ~ 6400 . If its dynamics were 
fully equivalent to the isothermal dynamics of a system of 
age i^eff ~ 6400, the PDFs of magnetization differences 
SM calculated over intervals of the same length St would 
overlap in the two cases. In Fig. ([5]), the PDF of isother- 
mal fluctuations for a system with t w = 6400 (stars), are 
compared to the corresponding PDF of a T-shifted sys- 
tem of age t w = 1000 (circles). Both data sets are plot- 
ted on a logarithmic vertical scale. The Gaussian parts 
overlap, but the intermittent tails do not. They are how- 
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FIG. 3: (Color on line) For the temperatures indicated, the four panels depict scaling plots of the average magnetization per 
spin, shifted by the constant bo(t w ). In each panel, three sets of data are shown. The ordinate is the average linear response 
under different conditions. The abscissa is, a part from an additive constant, the system age t, scaled with either the true or 
the effective age. The two nearly overlapping data sets are for: i) Isothermal response at the indicated temperature, versus the 
age scaled by t„. ii) T-shifted response, plotted versus the system age scaled by t w ,eS ■ The quality of the collapse of these 
two data sets gauges the relevance of the effective age parameterization, Hi) Same data as ii), but using the actual t w value to 
scale the age. 



ever parallel, meaning that the weight of the intermittent 
fluctuations relative to the weight of the Gaussian fluc- 
tuations is higher in the T-shifted than in the isothermal 
case, even though the normalized size distribution of the 
intermittent fluctuations is the same in both cases. Since 
the Gaussian fluctuations have zero average, the average 
response is described via £ W)e ff as previously discussed. 

In a record dynamics scenario, the marginally stable 
attractor reached at t w is associated to a dynamical bar- 
rier of magnitude b(t w ), which matches the largest among 
the thermal fluctuations experienced by the system dur- 
ing the aging process. Thermal fluctuations are described 
by the equilibrium Boltzmann distribution P(E) charac- 



terizing a single thermalized domain. In our case, the 
local density of states can be treated as a constant, and 
P{E) ex exp(—E/T). The typical value of the largest 
among t w fluctuations at temperature T then scales as 
b(t w ) oc Tlog(fw) WA - Note that the Arrhenius es- 
cape time corresponding to b(t w ) equals t w . j48|. Af- 
ter a negative T-shift from Ti to Tt < Ti, the Arrhe- 
nius time corresponding to the barrier b(t w ) becomes 

tw, e ff = exp(6(£ w )/T/) = tw^ Tf ■ We thus arrive at the 
conclusion that 

t w , e s = £w (T " T/) ; x(T h T f ) = (4) 

J f 
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FIG. 4: (Color on line) The four panels depict, for the temperatures indicated, scaling plots of the average magnetization per 
spin, vertically adjusted by the constant bo(t w ). In each panel, three sets of data are shown. The ordinate is the average linear 
response under different conditions. The abscissa is, a part from an additive constant, the system age t, scaled with either the 
true or the effective age. The two nearly overlapping data sets are i) isothermal response at the given temperature, versus the 
age scaled by t w . ii) T shifted response. The data are plotted versus the age scaled by t w ,eg. The quality of the collapse of 
the first two data sets gauges the relevance of the effective age parameterization. Hi) Same data as ii), but using the actual t w 
value to scale the age. 



Note that the temperature dependence of the exponent 
x given by the above formula has no adjustable param- 
eters. In Fig. ([S]), left panel, the empirically determined 
effective age is plotted versus the actual age on a log- 
log scale, for T-shifts 1.6 — > 1.5 (red squares) 1.8 — > 1.5 
(green polygons) and 1.6 — > 1.4 (blue diamonds). The 
corresponding lines arc obtained by linear regression, and 



their slopes provide a numerical estimate of the expo- 
nent x. Several other T-shifts are also considered, but 
the corresponding data are not shown for graphical rea- 
sons. The values of t%, corresponding to all available x 
values are plotted versus the empirical effective age in 
the right panel of Fig. [S] The color and symbol coding 
is as follows: 1.6 — > 1.5 (green circles), 1.7 — > 1.5 (red 
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FIG. 5: (Color on line) The PDF of magnetic fluctuations 
for isothermal (stars) and T-shifted (circles) protocols. The 
data are collected in intervals of type [t, 50t]. The initial ages 
are i w ,cff = 6400 for isothermal aging and t w = 1000 for T- 
shifted aging. The 8M values are calculated over intervals of 
size 6t = 500. A full collapse would imply the full equivalence 
of the T-shifted response and the isothermal response of a 
system of the correct effective age. The lines are fits to a 
zero centered Gaussian and to two exponentials, which are 
obtained using different subsets of the data. 

squares), 1.8 — > 1.5 (blue diamonds) 1.9 — > 1.5 (black 
hexagons), 1.6 — > 1.4 (green polygons), 1.8 — > 1.6 (red 
stars), 1.8 — * 1.6 (blue plusses) and 2.0 — > 1.8, (black 
upper triangles). The full line is the theoretical predic- 
tion given in Eq. ([4J. Thus, the distance between the 
data points and the line represents the mismatch between 
data and theory. Except for shifts involving a relatively 
high initial temperature, especially for high values of the 
effective age, the agreement is good. 

IV. DISCUSSION 

Summarizing the results of Refs. [H, [3(| EO, EH, us- 
ing the number of quakes n as an effective 'time' variable 
produces a translationally invariant rc-paramctcrizcd ag- 
ing description. This description relies on the statistical 
independence of the quakes and on the statistical subor- 
dination of physical changes to the quakes. It does not 
rely on the amount of change induces by the events. The 
(stochastic) changes may depend on the system type, on 
the observable considered, and even on n itself. They 
can either be dealt with with further modeling assump- 
tions EEl , or using the generic property that physical 
averages admit eigenvalue expansions, when treated as a 
function of n [3fl |. Averaging such expansions over the 
Poisson distribution of n turns each exponential term 
e.g. ae cn , into a fractional power (t/t w ) Xc , where the 



exponent A c is related to the eigenvalue c of the time 
re-parameterized problem by 

X c (T)=a(N)(e< N ^ -1). (5) 

The dependence of the exponent on the system size N 
and temperature T is now explicitely introduced. Note 
that the symbol A c here denotes a generic exponent, and 
is not to be confused with the (nearly) homonymous dy- 
namical exponent A widely used in the literature. Consid- 
ering that a{N) cx N, several exponents may exist which 
diverge with N. These exponent are dynamically irrele- 
vant as the corresponding modes quickly decay to zero. 
Some exponents may remain bounded in the thermody- 
namical limit, namely if the eigenvalue of the time- re- 
parameterized problem, c(N) vanishes as 1/N or faster, 
as N — » oo. Eigenvalues approaching zero faster than 
1/N correspond to frozen modes, and are also dynam- 
ically irrelevant. Observable power-laws are thus con- 
nected to relaxation 'time' scales of the re-parameterized 
dynamics which diverge linearly with system size, see e.g. 
Ref. [45| for an explicit demonstration of such linear di- 
vergence. The corresponding exponents have the form 
A C (T) = a(N)c(N,T), which is asymptotically indepen- 
dent of N for large N. 

Far away from equilibrium, and hence far from satura- 
tion, the smallest relevant exponent can be small enough 
to justify the further expansion (t/t w ) Xa oc 1 + A C ln(t/t w ) 
for a wide range of t/t w . If a logarithmic time dependence 
exists, it will be dominant for large values of t/t w . This 
is widely observed and implies that linear response and 
autocorrelation are again proportional. An effective tem- 
perature can then be defined, albeit in a non-universal 
manner [28| via the Fluctuation Dissipation ratio.. 

Numerical studies of domain growth convincingly show 
that the characteristic linear size of a thermalized domain 
increases algebraically in time [l,[l8j]. In this context, the 
T dependence of the growth exponent A indicates that 
the energy barrier surmounted at age t grows logarithmi- 
cally with t. The proposed explanation is however ad hoc 
and differs from the original assumptions of Fisher and 
Huse [47|, who expected a power- law scaling. In record 
dynamics, the same logarithmic scaling flows from the 
proportionality between said barrier and the typical size 
of extremal thermal fluctuations at age t. This property 
leads directly to the temperature dependence of the ex- 
ponent x given in Eq. Thus, even though the linear 
domain size R(t) does not appear conspicuously in the 
present description, central aging features are clearly as- 
sociated to properties of localized domains. Firstly, the 
presence of non-interacting domains can be deduced from 
the statistical properties of energy and linear response 
fluctuations [23, [U, [HJ , e.g. it produces the linear sys- 
tem size dependence of the factor a. Secondly, an approx- 
imately hierarchical energy landscape is required: In such 
landscape energy fluctuations of record size are, by con- 
struction, necessary to trigger any changes of metastable 
attractor. Last but not least, the temperature indepen- 
dence of a" [22|, HH is only possible if the landscape is 
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w w, eff 



FIG. 6: (Color on line) The left panel shows the algebraic relationship between the effective and the actual age for the three 
different T-shifts 1.6 — > 1.5 (red squares), 1.8 — > 1.5 (green polygons) and 1.6 — > 1.4 (blue diamonds). Other cases analyzed 
are not shown for graphical reasons The full lines correspond to the analytical prediction t w ,eff = £w, where x is determined 
by a least square fit. In the right panel, t%, is plotted versus £ w ,cff using different symbols and colors for all available x values. 
Deviations from the line gauge the discrepancy between predicted and observed behavior. 



invariant to a change of energy scale. 

In conclusion, a record dynamics description accounts 
for: i) the linear response phenomenology of aging sys- 
tems, including the effect of temperature shifts, ii) The 
origin of asymptotic descriptions in term of effective tem- 
peratures, in) The interplay of real space properties, i.e. 
domains, and hierarchical properties associated to the 
energy landscape of each domain. 
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